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ABSTRACT 

We present a series of simulations of the self-regulated growth of supermassive black holes (SMBHs) 
in galaxies via three different fueling mechanisms: major mergers, minor mergers, and disk instabil- 
ities. The SMBHs in all three scenarios follow the same black hole fundamental plane (BHFP) and 
correlation with bulge binding energy seen in simulations of major mergers, and observed locally. 
Furthermore, provided that the total gas supply is significantly larger than the mass of the SMBH, its 
limiting mass is not influenced by the amount of gas available or the efficiency of black hole growth. 
This supports the assertion that SMBHs accrete until they reach a critical mass at which feedback 
is sufficient to unbind the gas locally, terminating the inflow and stalling further growth. At the 
same time, while minor and major mergers follow the same projected correlations (e.g., the Mbh — c 
and Magorrian relations), SMBHs grown via disk instabilities do not, owing to structural differences 
between the host bulges. This finding is supported by recent observations of SMBHs in pseudobulges 
and bulges in barred systems, as compared to those hosted by classical bulges. Taken together, this 
provides support for the BHFP and binding energy correlations as being more "fundamental" than 
other proposed correlations in that they reflect the physical mechanism driving the co-evolution of 
SMBHs and spheroids. 

Subject headings: galaxies: formation - galaxies: evolution - black hole physics - galaxies: general - 
methods: numerical 



1. INTRODUCTION 

Over the past two decades, it has been established 
that supermassive black holes (SMBHs) in the nuclei 
of galaxies are common, and that their masses are 
connected to the properties of their hosts. As the 
sample of robust SMBH measurements has in creased 
(jKormendv fc R ichstond Il995t iKormendv * '2004") several 
black hole - host galaxy correlations (see Nov ak et al.l 
l2006l for a review) have been proposed, includ- 
ing c entral velocity dispers io n (iFerrarese fc Merritlj 
120001: iGebhardt etal] l2000al: iTremaine et all I2002D. 
bulge mass /luminos ity (jKormendv & Richstond 
19951 iMagorrian ct al. Il998t iMcLure fc Dunlop 200l 



process, along with gas dissipation, n aturally account 
for t h e observed SMBH corre la tions (IDi Matteo et al.l 
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concentration ( Graham et al .1 120011: iGraham fc Driveil 
|2007( ). and binding energy (jAller fc Richstond l2007t ). 
Recently, it has been suggested that these observed 
correlations may be projecti ons of a "B l ack-Ho le 
Fundamental Plane" (BHFP: [Hopkins et al.l I2q07allbl : 
lAller fc Richstond l2007t iBarwav fc Kembhavil \20m . 
analogous to that describing the structural prop- 
erties of elliptical ga laxies (jPressler et al.l 119871 : 
iDiorgovski fc Daviall987[ ). that is manifest as a tilted 
bulge binding energy correlation. These relationships 
are indicative of an intimate connection between SMBH 
growth and galaxy formation/evolution. 

There is also increasing evidence, both direct and 
indirect, that the growth of SMBHs is self-regulated, 
with strong feedback from the active galactic nucleus 
(AGN) playing an important role. Numerical simulations 
of major mergers of gas-rich spirals that include this 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138 
^ jyounger@cfa.harvard.edu 
^ Keck Foundation Fellow 



liptica l galaxies (Robertson ct al. 2006b; Hopkins et al.l 
|2008h), the struct ure (IHo pkins ct al. 2008c,f^c) and 
stellar kinematics ()Cox et al.. .20063 of massive ellip- 
ticals, and the redshift evo lution of the correlations 
(jHopkins et al.ll2007al [2008gi) . Moreover, a strong feed- 
back mode has proven essential in semi-analytic mod- 
els and numeri cal simulations of hierarchic al galaxy 
formation (e.g.. ICroton et al.l 120061: iBower et al.. ,200(1 
iSiiacki et all 120071: IDi Matteo et all |2008|, Sommerville 
et al. 2008, in prep.) and the merger- driven 
evolut ion of quasa rs and r ed galaxies (Hopkins et al.l 
[2005allbl lclldl. [2006al l3. [2008alldh . The prevalence of this 
process is further supported by direct observational 
evidence for outflows in both S eyfert galaxies (e.g., 
Crenshaw et al.' '1999^; 'Kriss et al.' '20001: iKaastra et all 
2000, 2002; Kaspi ct al. 2000a, 2002) and more luminous 
quasars (e.g..lWevmann et al.lll991l : iKorista et al.lll993l : 
lArav et al.l l2()0l[) driven bv strong thermal and radia- 
tive feedback origina ting from the AGN (see review by 
ICrenshaw et al.ll2003D . 

Despite these observational and theoretical advances, 
the fundamental character of these SMBH correlations 
remains poorly understood. It appears clear the SMBHs 
are closely linked with the st ructural properties of their 
host galaxy's bulge (see e.g.. iNovak et al.ll2006l ). How- 
ever, an important question has not been adequately ad- 
dressed: to the extent that the growth of SMBHs is self- 
regulated, what property - or properties - of this bulge 
does the SMBH "see?" The more "fundamental" the 
scaling relation, the more it reflects the physical mech- 
anism driving; the co-evolution of SMBHs and bulges. 
ins et al.l (|2007a|) suggest that the final mass of the 
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Fig. 1. — Top: Relative separation of the two SMBHs o ver t he course of an Mp/Ms = 8 (left) and 4 (right) merger simulation. Included 
are results using both the accretion drag (black line; see ^ 12.31 for details) and repositioning (red line; see Johansson et al. 2008) methods. 
The circles indicate the time of the merger. BOTTOM: Same labeling for the relative velocities of the two SMBHs over the course of the 
simulation. 



SMBH is set by the depth of the local potential; the 
SMBH grows until feedback unbinds the local gas supply, 
abruptly terminating its growth. However, this hypoth- 
esis, while promising, has not yet been systematically 
tested. 

In this paper, we use hydrodynamical simulations to 
examine the self-regulated growth of SMBHs in three dif- 
ferent fueling modes: major mergers, minor mergers, and 
disk instabilities. Our aim is to investigate the physical 
mechanism that determines the final mass of the SMBH. 
It is organized as follows: in §[2] we review our methodol- 
ogy, in §[3] we describe the simulations, in §[3] we present 
our results, and in §[5] we discuss implications. 

2. METHODS 

2.1. Hydrodynamics 

In what follows, we describe a series of simulations 
run using PGadget2 (^pringel 2005), a massively par- 
allel smoothed pa£ticle_Ji^ (SPH: see 
also iLucvl 119771 : iGingold fc MonaghanI I1977D using the 
entrop y conserving formulation of ISpringel &: Hernguist 
(l2002r) . coupled with a Tree algorithm ( Barnes fc Hut 
Il986f ) to com pute the gravitational for ces between parti- 
cles (see also iHernguist fc Katj|1989t). We also i nclud e 
the effects of radiative cooling (see iKatz et all 119961 ). 
and employ a sub-resolution model for the interstellar 
medium (ISM) to describe star formation and super- 
nova feedback, tuned to reproduce the local Schmidt Law 
(jSchmidtl 119591 : iKennicuttI [1998), that pressurizes the 
star- forming gas via an effective equation of state qeos 
that interpolates between an isothermal {qeos — 0) and 



full m ultiphase ISM {qeos — 1; ISpringel fc Hernguisli 
l2003f ). This sub-resolution model is essential for simula- 
tions of gas rich systems, because it allows us to evolve 
di sks of arb i trary gas content with out fragmenting owin; 
to iToomre (I1964D instabihty (see Springel fc Hernguis 
I20031I2005I: TRobertson et al.ll2004 l2006al) . 



2.2. SMBH Growth and Feedback 

SMBHs in our simulations are represented by "sink" 
particles, which increase their mass at the rate Mbh 
via Eddington-limited Bondi-Hoyle-Lyttleton a c cretion 
movie fc Lxttletonl [T94lt IBondi fc Hovli [1941 IBondil 
[1952) of the surrounding gas. As the numerical prefactor 
for the accretion rate - Bondi accretion is a dimensionless 
scaling - we adopt an efficiency of a — 100 for all the sim- 
ulations used in this w ork (see also ISpringel et aLll2005[ 
I Johansson et al.ll2008f ). However, this choice matters lit- 
tle as the SMBH's growth is Eddington limited during all 
the times of interest. The bolometric luminosity of the 
SMBH is assumed to be Ltoi = eMc^, where e = 0.1 is 
the canonical radiative efficiency for thin-disk accretion. 
Feedback is implemented by coupling a fixed fraction of 
this luminosity (feedback efficiency r/; typically rj « 5%) 
to its environment by injecting thermal energy into the 
surrounding gas, weighted by the SPH smoothing ker- 
nel. The net ratio between the dimensional accretion 
rate (Bondi or Eddington) and the thermal feedback in- 
put into the gas is in some sense the only free parameter 
-a, e, and r/ are degenerate - in our feedback model, 
and that ratio is chosen to repro duce the normalization 
of the local Mbh - (t relation ijDi Matteo et all 120051 : 
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Fig. 2. — Growth of the SMBH ov er th e course of an Mp/Ms = 8 (left) and 4 (rig ht) merger simulation. Included are results using both 
the accretion drag (black line; see 12.31 for details) and repositioning (red line; see [Johansson et al.ll2008l) methods. The circles indicate 
the time of the merger. There is no systematic trend with the choice of either the accretion drag or repositioning methods. 

[Robertson et al.ir2006cl ): the slopes and relative normal- 
izations of the resulting correlations of SMBH mass with 
galaxy properties are generic, and are generally not tun- 
able. 

It may be argued that this specific method of im- 
plementing feedback - in this case an isotropic ther- 
mal coupling - is significant for our results, especially 
in the absence of constraints from a well-understood, 
detailed model for AGN feedback. It is worth noting 
that alternate feedback prescriptions hav e led to simi- 
lar re sults as those presented here (e.g., iMurrav et al.l 
|2005| ). However, more generally our isotropic thermal 
coupling to the surrounding gas is a good approxima- 
tion to a variety of feedback mechanisms, as it leads 
to a shock front which isotropizes and becomes well- 
mixed over physical scales smaller than those relevant to 
our simulations, and in tim escales smaller than the dy- 
namical time of the g alaxy (jHopkins fc Hernauistll200a 
iHopkins et"aLll2006bh .'^ The spatial (i.e., SPH smooth- 
ing length) resolution for the simulations analyzed here 
is 30 — 50 pc, and the dynamical time over those scales 



is tdyn(30 — 50pc) ~ 10' years. Numerical experi- 
ments have shown that shock fronts, even if initially 
highly beamed, will become spherical on an isotropiza- 

tion timescale of tiso ~ 6xlO'^(£'5i/p24)^^'^ years-where 
-E51 is the input energy in units of 10^^ ergs and p2A is the 
density of the ambient IS M in units of 10~^"^ g cm~'^ - 
and size scales of <, 10 pc (|Aval fc Piranll2001[ ). both of 
which are below scales probed by our simulations. There- 
fore, for the simulations presented here, we believe that 
an isotropic thermal coupling prescription for feedback 
offers a good approx imation (for further discussion, see 
ins et al.ll2006af ): an assertion that we will test in 
detail in due course. 



2.3. Black Hole Mergers 

In simulations of interacting systems the implemen- 
tation of SMBH mergers can potentially introduce a 
systematic bias in the final SMBH mass. Absent res- 

* This approximation is n ot valid in the case of "radio-mode" 
feedback l lCroton et al.|[200g) . in which MHD effects can maintain 
a coUimated jet on scale comparable to the galaxy. However, obser- 
vations suggest that black holes gain most of their mass in bright 
quasar modes ([Spltan 1982), and < 10% of these objects are radio- 
loud. Therefore, we do not believe that this will be the relevant 
feedback mechanism for the growth mechanisms presented here. 



olution constraints, we would expect the BH binary 
to e ventually harden via a combination of dy nami- 
cal (jMakino fc Funatol [200l iBerczik et al.l I2OO60 and 
hvdrodvnamic (^Es cala et al. l '2004') processes, at which 
point gravitational wave emission dominates the energy 
loss and the two SMBHs in -spiral and coalesce (e.g., 
iMilosavlievic fc Merrittl[2"001l ). However, because com- 
putational limitations prevent following this evolution in 
a galaxy-scale simulation - and gravitational radiation 
is not modeled at all - th e standard prescription adopted 
bv lSpringel eTall (pOOl merges the two SMBHs if they 
are within a smoothing length with a relative velocity 
less th an or equal to th e local sound speed. 

Like [Johansson et al.l 1)2008 ') , we found that this sim- 
ple criterion was not in of itself adequate to ensure rapid 
merging of the SMBHs at the end of the merger - particu- 
larly for the unequal mass interactions. However, rather 
than employ their repositioning scheme - wherein the 
SMBH is repositioned at every timestep at the minimum 
of the local potential - w e introduced a term accounting 
for Bondi accretion drag (|Edgaij|2004l ): 



F, 



drag 



MbhVo 



(1) 



where Voo is the initial infall speed of the gas in the Bondi 
approximation. This is a very a pproximate treatment 
(for a more detailed analysis, see iRuderman fc Spiegel 
I1971D . but captures the essential physics at a level that 
is adequate for our purposes. 

To investigate the sensitivity of our results to this 
choice of numerical implementation, we ran two repre- 
sentative simulations - an Mp/Ms — 8 and 4 interaction 
with initial gas fractions of fg = 0.8^ and 0.4 respectively 
and identical orbital parameters (see § 13.11 for details) - 
with both accretion drag and repositioning. In Figure [T] 
wc show the relative separations and velocities of the two 
SMBHs, and in Figure [2] we show the SMBH growth his- 
tory. There is no apparent systematic trend based on the 
particularly choice of implementation. 

Based on this simple experiment, we do not see a sys- 
tematic trend in the use of accretion drag over reposi- 
tioning. Rather, the choice of one or the other method 
introduces a scatter of ^ 0.1—0.2 dex, which is noticeably 

^ This high gas fraction was chosen to ensure that the primary 
had a significant gas fraction — in this case fg K, 0.3 — at the time 
of the SMBH merger. 
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TABLE 1 

Progenitor Disk Structural Parameters 





V200 


M200 
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ho 
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[km s^-^] 
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Sb 


160 


95 


9 


4.1 


3.3 X 10^ 


2.0 X 10^ 


Sc 


130 


51 


10 


3.2 


1.8 X 10^ 


1.1 X 10^ 


Sd 


100 


23 


12 


2.2 


8.1 X 10* 


4.9 X 10* 


Im 


80 


12 


12 


1.6 


4.2 X 10* 


1.5 X 10* 



smaller than both the simulation-to-simulation scatter, 
and even the uncertainty associated with the stochastic 
nature of SMBH accretion. Therefore, while we note that 
a more systematic study of SMBH merger prescriptions 
for hydrodynamical simulations is warranted, our choice 
of the accretion drag method does not appear to bias our 
results. However, we do prefer it over repositioning for 
its numerical stability; in some - admittedly rare - cases, 
and particularly when the mass ratio of the interaction 
is high and one galaxy dominates the local potential, 
the central SMBH in the secondary can decouple from 
its host galaxy and follows the gravitational field lines 
of the primary, resulting in a spurious SMBH merger. 
However, in general, seeing as the choice of method does 
not qualitatively change our results, and there is no clear 
physical argument to favor either, we consider it one of 
the systematic uncertainties in our modeling. 

2.4. Progenitor Disks 

The progenitor disk m odels are constructed follow- 
ing [SEringeretlan (|2005f l. Exponential disks of gas and 
stars - and optionally a compa ct bulge - are em bed- 
ded in a dark matter halo with a 'Hernguistl (|199Cl[ ) den- 
sity profile, as motivated by cosmological N-body sim- 
ulations (e.g.. iNavarro et al.l ll996t iBullock et alTl2001t 
iBusha et al.l I2005D . Their total mass is VioQ/{10GHo), 
where V200 is the maximum circular velocity at an over- 
density of 200pc, with a disk mass fraction of m^, disk 
gas fraction of fg — mg/(md + rng), and a bulge mass 
fraction of rub- The disk scale-length is computed as- 
suming a spin parameter of A = 0.033 — 0.05, again 
motivated by cosrnological N-body siniulations (e.g.. 
Cole & Lacey '1996 yyitvitska et al.l l2002t iMaccio et atl 
2007; Bett et al. 2 00711 

For this work, we use three sets of simulations: (1) 
major mergers, (2) a mass ratio series, and (3) isolated, 
bar-unstable disks (hereafter referred to as simply "un- 
stable disks" ) . The structural parameters of the progen- 
itor disks u sed in the major merger series are discussed 
in detail bv lRobertson et all (j2006d lbl). and we refer the 
reader to those references for the specifics. The initial 
disks used in both the mass ratio and unstable disk se- 
ries are descr ibed in Table [II and are similar to those 
used iu iYounger et al.l (1200711 

3. SIMULATIONS 
3.1. Mergers 

The interactions in the major mergers series, and in 
particular thei r orbital parameter s, are discu ssed in de- 
tail by Robcrt sonet al.l ( 2006b, c). ICox et all (2006), and 
[Hopkins et all jlOOTi); we refer the reader to these pa- 
pers for a detailed discussion. Each galaxy typically con- 
tains 6 X 10^ halo particles (0.2ft,^^ kpc softening), 2 x 10* 




M /M -2 

M„/M =4 



Time [G/r] 

Fig. 3. — SMBH growth histories for four different interactions 
with mass ratios of: Mp/Ms ~ 8 (dash-double-dot), 4 (dash— dot), 
2 (dash), and 1 (major merger; solid). The orbital inclination for 
each was i = 30° and the gas fractions of both initial disks 40% 
if 9 — 0-4). The major mergers clearly grow much more massive 
SMBHs than do the minor {MpMs ~4-8) mergers. 

bulge particles (O.l/i^^ kpc softening), and 4 x 10'' stel- 
lar disk particles (0.1/i~^ kpc softening). The two pro- 
genitor galaxies are placed on a zero-energy parabolic 
orbit, as motived by cosmological N-body si mulations 
(e.g.. lBensonl[2005l: Khochfar fc BurkeTtI [200I . We in- 
clude a range of gas fractions, from relatively gas poor 
{fg — 0.05) to pure gas disks {fg — 1), covering « 2.5 
dex in total mass - e xcluded the scaled disks from 
iRobertson et al.l ()2006cf ). which are intended to be repre- 
sentative of high redshift systems. We have furthermore 
varied the resolution by up to a factor of 128 as many 
particles to confirm that our results are robust. 

The mass ratio series (summarized in Table [2]) includes 
all permutations of the progenitor disk models in Table[TJ 
which results in primary to secondary mass ratios rang- 
ing from Mp/Ms « 8 to major mergers {Mp/Mg = 1). 
Both initial disks are dynamically stable over a Hubble 
time, and have disk mass fractions of rud = 0.05. The 
mass resolution of the halo component was ~ 4 x IO^Mq 
per particle {0.\bh~^ kpc softening), while the mass res- 
olution of the baryonic component (stars and gas) was 
~ 3.5 X IO^Mq per particle (0.05/i~^ kpc softening) - the 
particle counts are summarized in Table [21 To simplify 
the interpretation of the results, none of these models 
initially had bulges. As before, the two disks are ini- 
tialized on zero-energy parabolic orbits, with an impact 
parameter of ~ /iz?,pj where h^ p is the scale length of 
the primary disk. We also includ e five different orbital 
inclinations, as in lYounger et aH (|2007f ): i = 0° (pro- 
grade coplanar), 30°, 90° (polar), 150°, and 180° (ret- 
rograde coplanar). We consider three different choices 
for the gas fractions of the init ial disks: (1) a stan- 
dard gas-rich series (as in IRobertson et al1l2006cf ) where 
both initial disks have /„ = 0.4 over all mass ratios 
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TABLE 2 

Orbital and Initial Disk Parameters for the Mass Ratio Series 
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and inclinations, (2) a series with more gas-rich lower- 
mass (Sd and Im) secondaries (as suggested by observa- 
tions; see [Roberts &z Haynes 1994; Bell & de Jong 2003; 
ISchombert eirani200ll : iGeha et aLil2006( ) with fg^p = 0.4 
and fg^s = 0.8 for a prograde {i = 30°) and retrogade 
(i = 150°) interaction, and (3) a high gas fraction series 
with fg = 0.8 for both initial disks (which m ay be rep- 
resentative of higher redshift proge nitors; see lErb et al.l 
[200l lF5rster Schreiber eFall [20061 ). over all mass ratios 
and again for a prograde and retrograde interaction. 

The process by which SMBHs grow and self-regulate 
in major mergers is described in de t ail by several 
authors (see e .g., ISpringel et al.l 120051 : iHopkins et al.l 
l2005al [dl. l2006af ). Broadly, tidal torques during close 
passages between disks induce bar instabilities in the 
stellar distribution. These stellar bars then drain the 
gas within a critical radius of its angular momentum 
([Hopkins et al.l[2008bf ) causing it to flow towards the cen- 
tral regio ns, which in turn fuels a nuclear starburst (see 
also Barnes fc Hernguis^ l 19911: i Mihos fc Hernauistl[l99l . 
119961: iBarnes fc Herngu ist 1996) and grows the SMBH 
nearly at an exponential rate. Eventually, when the ac- 
cretion rate reaches a critical value, feedback terminates 
the gas inflow and drives a galaxy scale "superwind" that 
shuts down both star formation and SMBH growth. At 
the same time, the interaction dynamically heats the stel- 
lar component, transforming it into a pressure-supported 
spher oid: a "red a r id dead" bulge-dominated galaxy (see 
also BarnesI 119921 : iBarnes fc Hernquist[ 119921 : iHernquistl 
[T992. 19931). 

When the mass ratio is sufhciently high {Mp/Mg ^ 4; 
a "minor merger"), the interaction does not entirely 
destroy the initial disk (for a detailed discussion for 
disk survival see iHopkins et al.l l2008bf ). but instead 
moves a smaller fraction of the stellar mass into a 
central bulge component (for a more detailed dis- 
cussion of the kinematics of disk heating by minor 



merg er s, see e.g. iQuinn fc Goodmai] 119861: iQuinn et al.l 
1993; 'Walker et al .l 119961 : iVelazauez fc White! [l999l : 
Kazantzidis et al. 2 0071 ). And, in general the nuclear 
starburst and SMBH accretion do not exhaust the avail- 
able gas supply, leaving a substantia l fraction in the disk 
component of the remnant (see also lCox et all [20081 for 
images and a detailed analysis). 

These clear differences motivate a look at the growth 
history of SMBHs in our simulations, and their sen- 
sitivity to the parameters of the interaction. A de- 
tailed investigation into the efficiency with which mergers 
grow b lack holes, and their luminous lifetimes (see e.g., 
ins ct al. 2005b, d, 2006a) for different mass ratios is 
a topic worthy of more detailed study, which we postpone 
to future work. However, as a qualitative illustration, we 
show the growth histories of SMBHs in mergers of differ- 
ent mass ratios and orbital inclinations. 

In Figure [3l we see that at fixed inclination and pri- 
mary galaxy mass, major mergers {Mp/Mg '^2 — 1) grow 
SMBHs that are nearly an order of magnitude more mas- 
sive than thos e produced in ininor me rgers {Mp/Mg ^ 4). 
As shown by IHopkins et al.l ([2007af l and in subsequent 
sections of this work, the SMBH traces the binding en- 
ergy of the bulge. Therefore, the efficiency of SMBH 
growth in interactions of differing mass ratio is likely a 
consequence of the efficiency of dynamical heating and 
bulge formation during those encounters. 

Qualitatively, in a minor merger , the firs t pas- 
sage induces a bar, ig as inflow ([Hernquistl 119891 : 
I Hernguist fc Mihoi Il995l ) and - if it is an interpen- 
etrating encounter - violently relaxing approximately 
Ms of the disk stars, which together form a bulge (see 
IHopkins et al.l [2008b( l^. It also tends to put the sec- 

® This suggests that as the impact parameter increases, the re- 
sultant bulge will tend more towards a pseudobulge. In fa ct, si nce 
no galaxy is truly isolated, the unstable disks outlined in ii l3.2l can 
be thought of as the limiting case of low mass satellites and halo 
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Fig. 4. — SMBH growth histories for two different mass ratio interactions - a major merger Mp/Ms = 1 (right panel) and minor merger 
Mp/Ms ~ 8 (left panel).— for 40% gas initial disks { fg = 0.4) over five different inclinations: i = 0° (coplanar prograde; dash), 30° (solid), 
30° (polar; dot), 150° (dash— dot), 180° (coplanar retrograde; dash-double-dot). With the exception of exactly coplanar interactions, which 
are somewhat pathological, the final SMBH mass is not particularly sensitive to the orbital geometry. 



ondary on a mostly radial orbit, and strip a signif- 
icant fraction of its mass (iBullock fc JohnstonI l2005t 
lAbadi et al.]|2006t iStewart et al.ll2007D prior to final co- 
alescence. A low angular momentum final orbit, coupled 
with a central mass concentration will tend to damp out 
further bar-driven gas i nflows (se e e.g., iShen fc SellwoodI 
2004; Hozumi & HcrnqjIiillQOi iBournaud et al.l 120051 : 
Debattista et al.i i2006i . and discussion in ij 13. 2p . This, 
in combination with less dynamical heating owing to 
the lower energetics of the interaction, suppresses fur- 
ther bulge growth, and with it the final burst of SMBH 
activity seen in major mergers. 

In Figure [4l we examine the effect of orbital inclina- 
tion - prograde coplanar (j = 0°) to retrograde coplanar 
{i = 180°) - on both minor (left panel) and major (right 
panel) mergers. In general, compared to retrograde in- 
teractions, prograde encounters produce stronger tidal 
responses in the disk that induce bars and drive the gas 
inflo ws that fuel bulg e formation and SMBH growth (see 
e.g.. lCox et al.ll200l . We see this effect in our simula- 
tions: while the effect is more dramatic for major merg- 
ers, for the limiting case of a coplanar interaction, retro- 
grade orbits result in a significantly less massive SMBHs 
than do prograde. However, when the spins of the two 
galaxies are not aligned, the resulting SMBH is largely 
independent of orbital inclination. 

These trends with mass fraction and orbital inclination 
qualitatively motivate the aim of this work: given a par- 
ticular encounter, what mechanism sets the final SMBH 
mass? 

3.2. Disk Instabilities 

The simulations in the unstable disk series evolve vari- 
ants of the same disk models in isolation, the parameters 
of which are summarized in Table [31 The mass reso- 
lution and softening lengths of the particles were kept 
the same as with the mass ratio series, and the parti- 
substructure seeding the bar instabilities via fiyby encounters. We 
have run a limited experiment, varying the impact parameter for 
a small subset of interactions and find evidence that the structural 
properties of the resultant bulges do tend towards pseudobulges for 
very large impact parameters (including the normalization offset in 
^Ibh ~ see ^ 14 3l and l5.3l l. However, a systematic investigation 
into the dependence of the structural properties of the bulge on 
the orbital parameters of the interaction is outside the scope of 
this paper, and thus we defer it to future work. 



cle counts are summarized in Table HI The disk mass 
fraction is increased while holding the total mass fixed, 
which seeds global instabilities by increasing the self- 
gravity of th e disk relative to its kinetic energy from 
rotat ion (see lOstriker fc Peeblei 119731 : lEfstathiou et al.l 
Il982f) . For three initial disk models (Sb, Sc, and Sd), we 
consider two disk mass fractions {md = 0.08 and 0.10), 
and three gas fractions {fg = 0.4,0.6, and 0.8). For a 
subset of these simulations, we consider three different 
seed black hole masses - Mbha = 10^, 10^, and 10® 
Mq - and two different equations of state ~ qsos = 0.5, 
1.0. 

Stellar bars are found generic ally in num erical simu- 
lations of isolated disk galaxies (iHohll [l97ll ) . and repre- 
sent a global instability in which rotatin g disks swing- 
amplify spiral density wav e perturbations I Toomrj|1981l : 
iBinnev fc Tremaind 119871 ). They have been shown 
empirically to develop when the self-gravity of the 
disk is comparable to i t s kinetic energy of r otation 
(jOstriker fc Peebles! [l973l : lEfstathiou et al.lll982f) . Most 
disk galaxies in th e local universe have bars (e.g . , 
lEskridge et"al] 120001: iMenendez-Delmestre et"all 120071) . 
and - despite conflicting observational results - it seems 
clear that a substantial fraction of h igh redshift z ^ \ 
disks also show bar features (e.g., iJogee et"al] 120041 : 
ISheth et al.ll2008f). These bars h a ve be en shown both 
theoreticallv (|AthanassouIal Il992l |2000|) and observa- 
tiona lly to drive gas inflows and nuclear starbursts 
( Jog ee et al.lll99l l2002l [20051: iPetitpas fc Wilsm3[200l . 
[20031 ). As a result, they represent a potentially promis- 
ing mechanism for fueling the growth of nuclear SMBHs, 
which ha s at times been in voked in semi-analytic mod- 
els (e.g., [Bower et aLll2006[ ) as their dominant mode of 
growth. 

Since t he ons et of the instabil ity is a co Uisionless pro- 
cess (see lTooin rc 1981; Binncv fc Tremain c 1987), it de- 
velops first in the stellar distribution. The stellar bar 
then drains the gas within a critical radius of its angu- 
lar mo mentum ([Barnes fc Hernauistill996l : lHopkins et al.l 
l2008bD . driving it inwards fueling a nuclear starburst; a 
process analogous to merger-driven gas inflows in which 
the bar is induced by the ti me-evolution of the tida l 
field during close passages (see lMihos fc Hernauistl[T99^ . 
This nuclear starburst concentrates some of the mass 
of the initial disk at its center, forming a pseudobulge 
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TABLE 3 

Initial Parameters for the Unstable Disk Series 





N&mc 


Disk Kdciss Frjiction 
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H H ,% / ^^^(S> ) 
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1 n 




U O^-jUO 
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BH Seed Mass 
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0.4 
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TABLE 4 

Resolution of the Unstable Disk Series 



Name 


Nhalo 


^baryons 


uSbOS 


3.2 X 10^ 


3.2 X 10^ 


uSc08 


1.7 X 10^ 


1.7 X 10^ 


uSdOS 


7.9 X 10"* 


7.8 X 10* 


uSbIO 


3.2 X 10^ 


4.0 X 10^ 


uSclO 


1.7 X 10^ 


2.1 X 10^ 


uSdIO 


7.7 X 10"* 


9.8 X 10* 




0.0 0.5 1.0 1.5 2.0 2.5 

Time [Gyr] 



Fig. 5. — Illustration of the bar-driven gas inflows and resulting 
feedback driven outflows in a representative unstable disk simula- 
tion: uSbOS with baryon fraction = 0.08 and initial gas fraction 
/g = 0.4. Shown is the SMBH mass (solid line) and total gas mass 
(Mg) within 100 (dashed Une), 250 (dotted line), 500 (dash-dotted 
line), and 1000 (dash-triple-dotted line) parsecs. The stellar bar 
efficiently concentrates a significant gas reservoir - several hun- 
dreds xMbh - within a few hundred parsecs of the SMBH which 
are subsequently expelled via feedback driven outflows. 

(see iKormendv fc Kennicuttl l2004al for a review) and 
eventually leading to the destruction of the bar over a 
relatively short timescale as the central mass concen- 
tration damps out the spiral density waves (~ 1 — 2 
Gyr; |Sh cn & Sellwoodi r2004: Hozumi'& Hcrn auistl 120051 : 
iBournaud et al.ll2005t iDebattista et al.ll2OO60 . Because 
there is no interpenetration by a satellite, there is corre- 
spondingly no violent relaxation of the stellar disk and 
as a result the pseudobulge is almost entirely dominated 
by stars formed during the simulation. 

Though the feedback model employed for the unstable 
disk series is identical in both implementation and pa- 
rameter choice to that used for the major mergers and 




Time [Gyr] 

Fig. 6. — Local gas density - the mean density over one SPH 
smoothing length of 64 particles — near one of the SMBHs in four 
merger simulations of varying mass ratio - Mp/Ms = 1 (black 
line), 2 (yellow line), 4 (green line), and 8 (red line), all with an 
orbital inclination of i = 30° and initial gas fraction fg = 0.4 — 
as compared to a representative unstable disk simulation (same as 
in Figure [5] blue line). The local gas densities, which are used 
to calculate the instantaneous accretion rate of the SMBH, are 
comparable in mergers and unstable disks. 

mass ratio series, in order to insure a fair comparison 
it is important that the growth of the SMBH be self- 
regulated by feedback rather than limited by gas de- 
pletion in the central regions. This will be the case if 
the stellar bars, like those induced during close passages 
during major and minor mergers, drive gas inflows that 
create a central gas reservoir many times more massive 
than the SMBH and a substantial density enhancement 
in its immediate vicinity, which will result in Eddington- 
limitcd growth until the local gas supply is then expelled 
by a thermal feedback driven blast wave. In Figure [5] 
we show the gas inflow - and feedback driven outflow - 
as a function of time for a representative unstable disk 
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simulation. We find that the gas inflows driven by the 
stellar bar concentrates a significant gas supply in the 
central several hundreds of parsecs. Furthermore, in Fig- 
ure [6] we confirm that the gas density enhancement in 
the immediate vicinity of the SMBH ~ which is used to 
calculate the instantaneous accretion rate ~ is compa- 
rable to those produced in merger-driven inflows. Fi- 
nally, we find that the accretion rate - i.e., AGN light 
curve - decays rapidly with a power-law slope similar 
to that expected f rom feedback driven outflows (e.g., 
iHopkins fc Hernouist 2006). As a result, we confirm that 
SMBH growth in our unstable disk simulations is regu- 
lated by feedback rather than gas depletion, and thus 
- if we assume that accretion and feedback in mergers 
and unstable disks are similar - our feedback model is 
appropriate to these systems. 

3.3. Systematic Uncertainties 

[Robertson et all (|2006cD and iHopkins et all (|2007aD 
both comment on the sensitivity of the predicted SMBH 
scalings arising from major mergers to choices in our 
models: in particular the choice of the initial seed mass 
Mbh,! and the EOS parameter qeos, which unlike fg do 
not have a clear observationally preferred, physically mo- 
tivated value. These authors find that over a relatively 
wide range of values for these parameters, the scaling 
relations with final SMBH are well within the scatter 
both in the observed and simulated relations. Therefore, 
while the overall growth history of an individual simu- 
lation may change significantly with specific choices for 
these parameters, the predicted scalings remain statisti- 
cally unchanged. This is consistent with our interpreta- 
tion of the scaling of SMBH mass with the binding energy 
of the bulge as the more fundamental, physical relation; 
the parameters of the bulge component of the remnant 
are relatively insensitive to these parameters. 

Because the mass ratio series is consistent with the 
scalings derived from major mergers, we are confident 
that they too are not affected by these systematic un- 
certainties in our modeling. However, it is not clear a 
priori how varying MsH^i and Qeos will affect SMBHs 
grown in unstable disks. This is because the SMBH does 
not grow by more than one or two orders of magnitude 
and the stochastic nature of its accretion may make it 
more sensitive to the details of the growth history than 
SMBHs grown via major or even minor mergers, which 
accrete several orders of magnitude more at comparable 
total baryonic mass (see, e.g. Figure [3]). 

In Figure [7] we show the growth history of SMBHs 
in three different disk models with fg — 0.4 for dif- 
ferent seed masses: uSbIO {Mbh,z = 10^ IO'^Mq/Zi), 
uSclO iMBH,r = 105,106Mo//i), and uSdIO (Msff,» = 
lO'^, lO^Mo/ft.'). We find that, as with the gas supply, 
while the detailed evolution differs substantially, the final 
SMBH masses are largely consistent to within the scatter 
in the simulated and observed relations. The largest dif- 
ference (~ 0.3 dex) is for the lowest mass disk model, for 
which the SMBH accretes the least mass and is therefore 
more sensitive to certain choices of parameters for our 
modeling. 

In Figure in we show the SMBH growth history for dif- 
ferent effective equation of state parameters qeos — 0.5 
and 1.0 for three different disk models with fg = 0.4. 
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Fig. 7. — SMBH growth histories for three different disk models 
- Sb (solid line), Sc (dashed line), and Sd (dot— dot— dashed line) - 
and three different SMBH seed masses — 

MBH,i = 10" (blue), 10^ 
(black), and 10^ (red) Mq. We find that, despite differences in 
the SMBH growth histories, there is no systematic trend between 
the final SMBH mass and Mgu i, and that the maximum deviation 
for an order of magnitude difference in the seed mass is only ~ 0.3 
dex. 




0.0 0,5 1.0 1.5 2.0 2,5 3,0 3,5 
Time [Gyr] 

Fig. 8. — SMBH growth histories for three different disk models 

- Sb (solid line), Sc (dashed line), and Sd (dot-dot-dashed line) 

- and three different SMBH seed masses — Qeos = 0.5 (red) and 
1.0 (black). We find - as with the initial seed SMBH mass - de- 
spite differences in the growth histories there is no systematic trend 
between the final SMBH mass and the stiffness of the EOS. 

As qEo s approaches unity, it incre a ses th e dynamical 
stab ihty of the gas agains t iToornrd (|1964f ) instabilities 
(see [Robertson et al.l^Op. l2006d). This effectively de- 
lays the onset of the bar-driven inflow in disk models 
with qEOS = 1-0 relative to those with a softer EOS 
Qeos = 0-5- However, despite these differences in the 
growth history of the SMBHs, their final masses are 
roughly consistent with one another to within the scatter 
in the observed and simulated relations. The largest dif- 
ferences, as with different seed masses, are for the lowest 
mass disk model 0.3 dex). More important, there are 
no systematic trends with qeos- the final SMBH mass 
in uSblO is somewhat lower, while those in uSclO and 
SdlO are shghtly higher. 

3.4. Remnant Properties 

The structural and kinematic properties of the rem- 
nant were extracted after the system reached a state of 
approximate dynamical equilibrium; for the major and 
mass ratio series, this was typically ^ I/i~^ Gyr after the 
final coalescence of the two SMBHs, while for the unsta- 
ble disks this was typically ~ 0.5h^^ Gyr after the peak 
of the SMBH accretion rate. If the remnant was disk- 
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Fig. 9. — Left: Binding energy correlation — which is statistically equivalent to the BHFP fH07: [Hopkins et al.ll2007al ) - for the major 
and mass ratio series (solid line, scatter indicated by dotted lines). The grey hexagons are the major mergers series and the grey inverted 
triangles are observations, both from H07. The mass ratio series — including all gas fractions - are colored according to the mass ratio (see 
color bar to the right of the figure) with the primary galaxy model indicated according to: Sb (filled circle), Sc (filled triangles), and Sd 
(filled square). Right: The binding energy correlation for mergers (see previous; the mass ratio series as filled triangles) and unstable disks 
(same labeling as primary galaxy in previous) for three different gas fractions: fg = 0.4 (red), 0.6 (green), and 0.8 (blue). We find that 
SMBHs grown via major mergers, minor mergers, and disk instabilities all lie along the same BHFP and binding energy correlations. 
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Fig. 10. — Influence of gas fraction on the location of SMBHs on 
the binding energy correlation for the mass ratio series. The label- 
ing is the same as the left panel of Figure [9] Shown are the mass 
ratio series simulations with initial gas fraction fg = 0.4, and the 
arrows indicate the location of identical interactions with more gas 
rich progenitors {fg = 0.8). We find that increasing the gas con- 
tent of the progenitors drives these systems along these relations 
- owing to the structural properties of bulges formed during more 
gas— rich interactions (see e.g.. [Robertson ct al. 2006a, b) - but not 
systematically away from them. This reflects the different struc- 
tural properties of bulges formed from more gas— rich disks, leaving 
the fundamental character of the BHFP and binding energy corre- 
lations unchanged. 



dominated, w e then fit a combination exponential and 
iSersid ()1968() profile to the projected stellar mass den- 
sity (i.e., "surface brightness") profile along 64 lines of 
sight^, from which we estimated the median bulge mass 
(Mb) and effective radius (Re). If the remnant was bulge- 
dominated (B/T ^ 0.8), we took Mi, to be the total stel- 
lar mass, and Re the half-mass radius. The inner velocity 
dispersion (a) was measured within Re, again taking the 
median along 64 different lines of sight. This procedure 
was identical between the mergers and unstable disks. 

''' The detailed formation and growth of stellar bulges via minor 
mergers is certainly an interesting problem in its own right, and 
the simulations presented here are ideal for just such an analysis. 
However, for this work we restrict ourselves to simple profile fits 
and bulge masses for both the purposes of clarity and to best ap- 
proximate observational estimates of bulge masses, and postpone a 
detailed analysis of the bulge structure and kinematics to a future 
paper (Cox et al. 2008, in preparation) 



Furthermore, iHopkins et al.l (|2008b[ ) show that surface 
brightness fitting of the kind employed in this work re- 
covers on average the bulge-to-disk ratios inferred from 
a kinematic decomposition of the stellar particles. 



4. RESULTS 



4.1. 



The Black Hole Fundamental Plane 

[Hopkins et al.l (|2007af ) find that there exist correla- 
tions in__tlie_j;esidualso£fits such as the Mbh — cr 
and iMagorrian et al.l (|1998( ) relations; for example, at 
fixed a the residual SMBH mass scales approximately as 



M°-^^, and at fixed M}, the residuals scale as 



^1.40 



By marginalizing over two parameters - Mt, and either 
(J 01 Re ~ they found a best-fit BHFP (analogous to 
that observed for elliptica l galaxies; [Dressier et al.|[l987l : 
[Diorgovski fc David [1987f ) which minimized these resid- 
ual correlati ons with the form Mb h ^ A/^-^^cr^'^, both in 
simulations (Hopkins et al.''2007a) and for observed sys- 
tems (jB arwav fc Kembhavi 2007; Hopkins et al. 2007bj; 
[Aller fc Richstond [20071. This is statistically indistin- 
guishable from a correlation with the bulge binding en- 
ergy proxy MfoCT^ (Mbh ~ (A^&cr^)° ''), and can be in- 
terpreted as reflecting the nature of feedback regulated 
SMBH growth: accretion accelerates until feedback is 
sufficient to unbind the local gas supply, abruptly termi- 
nating the inflow and cutting off further growth. There- 
fore it is more "fundamental" than its various projec - 
tions, such as the Mbh — f and IMagorrian et al.l ( 19981 ) 
relations - a point we discuss in more detail in § [5.1[ 

To systematically test this hypothesis, we examine the 
binding energy correlation - which is statistically equiva- 
lent to the BHFP - in three different modes of SMBH fu- 
eling: major mergers, minor mergers, and unstable disks 
(see § |3| for details) . Figure |5| shows the binding energy 
correlation for major mergers (grey hexagons), the mass 
ratio series (left: colored points; right: black triangles), 
and unstable disks (rig ht: colored poi n ts), alo ng with the 
observations listed in iHopkins et al.l ()2007bD . We find 
that over a range of baryonic masses, gas fractions, and 
orbital parameters they all lie along the same relation to 
within the scatter, and that all reproduce the observed 
correlation. 
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Fig. 11. — Correlation between the SMBH mass residual at fixed binding energy and the gas fraction at the peak of the starburst (/g,pcak) 
for mergers (left) and unstable disks (right), with the same labeling as Figure |9] We find that these residuals are uncorrelated in all three 
sets of simulations, indicating that the final SMBH mass is not determined by the available fuel supply. 




Fig. 12. — Mbh — o" relation for mergers (left) and unstable disks (right), with the same labeling as the right hand panel of Figure l9l 
The black solid and dashed lines indicate the best-fit relations and scatter from H07, and we include obs ervations of SMBHs hosted by 
pseudobulges — which are likely grown via d isk instabilities (see IKormendv h Kennicuttll2004al ) — from lHiil (H08; grey stars indicate data, 
dashed line indicates best— fit relation: I2008IV 



4.2. The Role of the Fuel Supply 

Inasmuch as their growth is terminated at a criti- 
cal accretion rate, the final masses of SMBHs should 
not be determined by the available fuel supply so long 
as the gas reservoir is much more massive than the 
SMBH. While it is the case that the final SMBH mass 
is strongly correlated with the total gas mass of the sys- 
tem, this reflects the structural properties of bulges in 
gas-rich merger remnants: owing the effects of dissi- 
pation, a more gas rich progenitor will lead to a more 
compact bulge at fixed total mass, and consequently a 
dee per central potent ial 4)c and v elocity dispersion a (see 
e.g.. lRobertson et al.l[2b 06b.c: H opkins et al.ll2007a[ ). As 
shown by iHopkins et al.. (,2007a)) for major mergers, at 
fixed potential the gas fraction has no effect on SMBH 
growth. 

One method of illustrating this is presented in Fig- 
ure 1101 which shows the BHFP and binding energy cor- 
relations for the mass ratio series (colored points) as com- 
pared to those derived from simulations of major mergers 
(grey hexagons; solid line with the dotted line indicating 
the scatter). Here we show that increasing the initial gas 
fraction from fg = 0.4 to 0.8 for the identical interac- 
tion will drive the remnant along the BHFP and binding 
energy correlations, but not systematically away from 
them. Therefore, we find that the typically more mas- 



sive SMBHs resulting from simulations of more gas-rich 
encounters is related to the structural properties of the 
resulting remna nt - which is more compact at fixed 
(|Hopkins et al.l i2008c.l.e.g,h). and thus has larger bind- 
ing energy - than it is to the larger available fuel supply, 
leaving the nature of these correlations unchanged. 

In Figure[Tll we show the residual correlation of SMBH 
mass at fixed binding energy with the gas fraction at the 
peak of the starburst (/g.poak) for both mergers (left) and 
unstable disks (right). We find no correlation between 
this residual and /g,peak- This demonstrates again that 
to the extent that SMBH growth is self-regulated, its 
final mass is not correlated with the available fuel supply 
but rather with the binding energy of its host bulge. 

4.3. The Projected Correlations 

While the BHFP and binding energy scalings with 
SMBH mass are consistent between mergers and unsta- 
ble disks, individual projections of these relations are not 
consistent. In Figures [T2l and fT3| we s how the predicted 
Mbh — tr and iMagorrian et al.l fl9981 scalings from the 
unstable disk simulations, as compared to both the major 
merger and mass ratio series. One of the more striking 
features of these results is the normalization offset in the 
Mbh — relation, in which SMBHs grown via disk in- 
stabilities are less massive at fixed a than those grown 
in mergers but maintain a similar slope. Or conversely. 
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Fig. 13. — IMagorrian et all |(T99 8) relation for mergers and un- 
stable disks, with the same labeling as Figure [T21 

the bulges grown via disk instabilities are more compact 

at fixed SMBH and bulge mass.^ 

This difference was noted by |h3 (|2008f) and iGrahamI 
((2008) as a normalization offset between the Mbh — c 
relations in "classical" bulges, which are likely the re- 
sult of mergers, and "pseudobulges" or bulges in barred 
systems which were likely grow n via secular processes 
()Kormendv fc Kennicuttll2004b[ ). Our simula tion s agree 
quite well with these observational results. iHul ()2008f ) 
interprets this result as possible evidence for less effi- 
cient fueling of SMBHs via bar instabilities, creating a 
less massive SMBH at fixed a. However, we find that 
this reffects structural differences between bulges pro- 
duced in merger simulations versus those produced via 
secular processes in the unstable disks: bulges and p seu- 
dobulges^ lie along different Faber & Jackson] l|1976( l re- 
lations (see Figure [H)) . These systems all lie along the 
same BHFP and binding energy correlations, suggest- 
ing that the lower mass SMBHs found in pseudobulges 
at fixed a reflect their lower overall binding energy ow- 
ing to structural differences rather than the efficiency of 
fueling. These structural differences are driven largely 
by the different formation mechanism of pseudobulges, 
which retain a great deal more rotation than their classi- 
cal analogues - the LOS velocity dispersion will include 
rotational motions. In addition, because all the mergers 
presented in this work coalesce, the structural differences 
can further be understood via simple energy arguments: 
in a merger, conservation of kinetic energy requires that 
the effective radius of the bulge increase - and therefore 
the velocity dispersion must decrease - when the system 
has fully merged (e.g.. iHernquist et al.lll993l ). 

5. DISCUSSION 

5.1. The "Fundamental" Character of the BHFP 

Among local systems, there are a relatively large 
number of tight correlations between SMBHs and the 
properties o f their host galaxies (for reviews, sec e.g., 
iNovak et all [2006: Hopkins ct al. 2007a; Graham 2008). 
The question of which of these is more "fundamental" 
has two components: (1) which is a better predictor of 
SMBH mass in systems for which detailed measurements 

* The ex act location of the pseud obulges produced in unstable 
disks on the lFaber fc JacksonI ||197(]|1 relation is somewhat sensitive 
to the initial conditions. However, their manifestly different struc- 
tural properties relative to classical bulges produced via mergers is 
a robust result. 



Fig. 14.— IFaber fc Jackson" (19761) relation for bul ges grown 
via mergers and disk instabilities, with the same labeling as the 
right hand panel of Figure |9] Bulges produced via major and 
minor mergers ("classical bulges") lie along the same correla- 
tion, wh ile those produced via secular processes ("pseudobulges"; 
IKormen dv fc Kcnnicutt 2004a) follow a different relation. This 
demonstrates that these two types of bulges are structurally differ- 
ent; plausible, given their likely different formation mechanisms. 

are not feasible, and (2) which better reffects the phys- 
ical mechanism driving the co-evolution of SMBHs and 
bulges. The answer to (1) is largely determined by obser- 
vational uncertainty and sample selection, and is outside 
the scope and limits of this investigation. However, (2) is 
directly addressed by our simulations, and highlights the 
physical mechanisms causing the self-regulated growth 
of SMBHs. 

IGrahamI (12001 argues that the residual-residual cor- 
relations that initially motivated the observed and sim- 
ulated BHFPs are a consequence of sample selection, 
which includes bulges produced via both interactions and 
secular processes. In particular, the author notes that 
when barred systems are excluded, these correlations are 
no longer significant. This is entirely consistent with ex- 
pectations from our modeling. Owing to the relatively 
small dynamic range probed by classical bulges with ro- 
bust SMBH mass measurements, we would not necessar- 
ily expect the residual correlations to appear significant 
without the inclusion of lower mass systems, which are 
more likely to be barred. More important, the correla- 
tion between residuals will be most apparent in systems 
with different structural properties, like classical bulges 
and massive ellipticals as compared to pseudobulges and 
barred systems. 

Moreover, ou r results suggest that while the results of 
IGrahamI ()2008[ ) directly address (1), they do not provide 
convincing evidence against the BHFP and binding en- 
ergy argument as the solution to (2). In Figure [H we 
show that simulations of three different fueling mecha- 
nisms - major mergers, minor mergers, and disk insta- 
bilities - all reproduce the observed BHFP and binding 
energy correlations. This suggests that in all three cases, 
the mass of the SMBH is set by the critical accretion rate 
at which feedback terminates the gas inffow, and halts 
further growth of the SMBH. At the same time, SMBHs 
grown in unstable disks versus interactions lie along dif- 
ferent projected correlations in Figures [12] and [131 In 
this way, the BHFP and binding energy correlations are 
more fundamental in that they reffect the feedback regu- 
lated growth of SMBHs across bulges with very different 
structural properties. 
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Fig. 15. — Left: Correlation between SMBH mass and the binding energy of all the stellar particles i?a,bindi for the mass ratio series. 
Labeling is the same as the left hand panel of Figure |9] The solid line represents a fit to the relation for major mergers {Mp/Ms < 2), 
which are bulge— dominated, with the r.m.s. scatter indicated by dotted lines. Right: The residual SMBH mass at fixed Ball, bind: 
compared to E^, ^^i^^/E^ii i,i„^ using the relation derived for major mergers (see left), where i?bulgG,bind ~ 10.1Mf,(T^ is the approximate 
total binding energy assuming a spherical I Hern quistI 111993 1 profile. We find a significant residual correlation as E^^ bind/^aii bind approaches 
unity, suggesting that the bulge binding energy is a better predictor of final SMBH mass than the total stellar tjinding energy in systems 
that retain a significant disk component (minor mergers; Mp/Ms <; 4); i.e., the SMBH "sees" the bulge more so than the overall stellar 
distribution. 

5.2. Implications for SMBH Fueling Models 

Our modeling demonstrates that in feedback-regulated 
models of SMBH growth, while the BHFP and bind- 
ing energy correlations are univers al, their projections 
- Mbh — iMagorrian et al.1 (|1998D . etc. - are sensitive 
to the structural properties of the spheroidal component 
of the host galaxy. For the particular case of Mbh — c, 
this results in a normalization (see § 14.31 and Figure [T^ 
offset between classical bulges - which were likely formed 
during mergers - and pseudobulges and those found in 
barred systems - which are likely produced via secular 
processes. Given that our particular implementation of- 
fers a good approximation to the effects of feedback more 
broadly (see § 12. 2|) , this is a generic prediction of models 
in which SMBH growth is self-regula ted via feedback. 

Co nsequently, that both th e BHFP (|Aller fc Richstond 

l2007t iHopkins et alJ l2007bl ) and an apparent offse t in 
Mbh — C7 between classical and secular bulges ()Hul 
l2008t lGraha^l2008h are observed provides evidence for 
feedback regulated inodels more generally (ISilk &: Rees 
1998l:lKine'20 0,l [20051: iMurrav et al.l2nO,ql:lSazonov et al " 



2005; Rober tson et all l2006cl: iThacker et al.l 120061 : 
Hopkins et al.n2007a[ ). It also presents a challenge 



to alternative models in which the projected scalings 
are relatively insensitive to the structural properties of 
the spher oid, including: s tellar capture by the accre- 
tion disk (jZhao et al.ll2002t iMiralda-Escude fc Kollmeierl 
I2005D regulation via a vis cous star-forming accretion 
disk (iBurkert fc Silk! |2001[). a diabatic SMBH growth 
(iMacMillan fc Hcnrikse^ [IqoI) , and direct gas collapse 
Adams et al..,200L ,20031 ). 



5.3. Implications for Evolution in the Mbh — cr Relation 

The existence of these two complementary Mbh — cr 
relations also has potential consequences for measure- 
ments of the evolution of this correlation. Studies of 
AGN at intermediate redshift have led to somewhat am- 
biguous conclusions: some show evolution in M bh ~ o" 
focu ct al. 2004, 2007; Wahcr et al. 2004; Wo o et all 
[200a 2008), while others do not (jShields et al.ll2003f ). 
Measurements of SMBH masses at intermediate red- 
shift are done using reverberation mapping techniques 



(jBlandford fc McKed [1981 iPetersonI [l99l . which use 
spectral variability to infer the size scale of the broad 
line region and - assuming circul ar orbits - estimate 
the mass of the S MBH (see also IWandel et all 119991: 
iKaspi et al.ll200dH l. These virial relations are calibrated 
using a relatively small sample of local , relatively low- 
luminositv AGN (iGebhardt et aLll2000b|: iFerrarese et all 
I200H lOnken et al.l l2004t iNelson et al.l I2004D . many of 
which may reside in pseudobulge hosts and therefore cal- 
ibrating them as classical bulges may not be appropriate. 

However, at cosmological distances, AGN will be pref- 
erentially more luminous, and therefore more likely to be 
fueled by mergers than disk instabilities. This could lead 
to inferred evolution in the Mbh — f relation which is 
more related to selection effects than true evolution in 
the correl ation. In fact, the shift of A loga ~ ~0.45 ob- 
serve d bv lTreu et al] (|2004 [2001 and I Woo et all (|2006l . 
l2008l ) is entirely consistent simply with the difference be- 
tween the normalizations of the merger and disk insta- 
bility Mbh — f relations in Figure [T^ 

5.4. Why the Bulge? 

Both observational and theoretical work has demon- 
strated that the SMBH is correlated with the structural 
properties of its host galaxy. However, these correlations 
are far more significant when applied to the spheroidal 
component - or bulge - rather than the galaxy as a whole. 
We illustrate this in Figure [151 in which we show the 
binding energy correlation including all the stellar parti- 
cles (ii'bind.an; left), and the residual SMBH mass at fixed 
-Ebind.aii as a function fo i^bind,aii/£'bind,b on the right, 
where £^bind,b is an estimate o f the total bulge binding 
energy assuming a iHernquis^ (|1990( ) profile. When we 
use the correlation for major mergers, there is a clear 
trend in the residuals - a ~ 5<t significant correlation - 
towards an overprediction of the final SMBH mass when 
the bulge mass is small; i?bind,an is not a good predictor of 
the SMBH for systems with significant disk contribution 
to the overall binding energy as compared to £'bind. bulge- 

The question then arises of why, even in minor merger 
remnants where it contributes a small fraction of the to- 
tal mass of the remnant, does the SMBH scale so tightly 
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Fig. 16. — Prom the major mergers series: a comparison of the 
scaling of the final SMBH mass with the binding energy of the 
stellar bulge particles after the remnant reaches a state of approx- 
imate dynamical equilibrium (black triangles) and gas particles at 
the peak of the starburst (grey triangles), along with power-law 
fits to each (solid and dashed lines respectively). We find that the 
slope of the gas binding energy is steeper than that of the stellar 
particles, indicating that the relationship between SMBH mass and 
bulge binding energy reflects the SMBH self-regulating its growth 
by terminating the gas inflow, and that gas forming the bulge stars 
that dominate the potential close to the SMBH. 

with the binding energy of the spheroidal component of 
the potential. There are two possible scenarios that may 
explain this correlation: (1) the bulge confines the gas 
inflow isotropically so that the SMBH can continue to 
accrete without ejecting gas along potential minima, or 
(2) the gas content of the final starburst dominates the 
local potential near the SMBH and therefore it grows un- 
til feedback is sufficient to unbind this local gas supply 
and terminate the inflow. While we defer a more thor- 
ough analysis of these two effects to future work, here we 
present some suggestive evidence for the latter over the 
former. 

In Figure [16] we present the SMBH scaling with two 
different binding energy measures: that of the stellar 
component (the "stellar binding energy" ) after the rem- 
nant has relaxed, and that of the gaseous component 
(the "gas binding energy" ) at the peak of the final star- 
burst. We furthermore restrict ourselves to the major 
merger series, for which it is a fair approximation that 
the remnant - including both the stars formed during the 
simulation and initialized with the progenitor disks - is 
entirely bulge dominated and largely gas free. We find 
that the relationship between the SMBH mass and the 



gaseous binding energy is somewhat steeper than that 
with the stellar binding energy: the logarithmic slopes 
are a = 0.85 ± 0.17 and 0.74 ± 0.13 respectively. More- 
over, this scaling with gas binding energy is statistically 
consistent with a linear relationship. This suggests - al- 
beit tentatively - that (2) is more likely the dominant 
process. However, since the slope is not quite linear, it 
also argues for some contribution from the confinement 
of gas by the spheroidal potential. 

6. CONCLUSIONS 

We present a series of simulations of the self-regulated 
growth of SMBHs in which Eddington-limited accretion 
continues until a critical accretion rate, at which point 
feedback is sufficient to terminate the gas inflow and cut 
off further growth. Simulations of three different fuel- 
ing mechanisms - minor mergers, major mergers, and 
disk instabilities - all follow the same B HFP and bind- 
ing energy correlations observ ed locally l|Hopkins et al.l 
I2007at lAller fc Richstond [20071 ). Increasing the gas con- 
tent of the initial disks shifts them along these correla- 
tions owing to structural changes in the resulting bulge 
arising from the effects of gas dissipation, but does not 
change their character. And, while the major and minor 
mergers both follow the same p r ojecte d correlations - the 
Mbh — cr and lMagorrian et al.l (|1998[ ) relations - the un- 
stable disks lie along different correlations, in agreement 
with observati ons o f pseudobulges and bulges in barred 
systems ((Hui .20081 iGrahamI [20081 ) . again reflecting the 
structural differences between b ulges formed via mergers 
and se cular processes (see also iKormendv fc KennicuttI 
l2004al) . Taken together, these simulations support the 
BHFP as the most "fundamental" scaling relation, in 
that it reflects the physical mechanism driving the co- 
evolution of SMBHs and bulges. 

Thanks in particular to the anonymous referee for their 
helpful comments on this manuscript. We also thank 
Marijn Franx, Dusan Keres, Stephanie Bush, Gurtina 
Besla, Chris Hayward, Nick Scoville, Dave Sanders, Lisa 
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work is supported in part by a grant from the W. M. 
Keck Foundation. 
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